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Summary 

Two computational methods are used to predict 
the flow over a generic helicopter fuselage of a simple 
configuration. A thin-layer Navier-Stokes code and 
a panel method code are used to compute the sur- 
face pressures for comparison with data from 4 ex- 
perimental conditions at 14 fuselage stations. The 
findings of both methods are in agreement with the 
experimental pressure data. However, separation 
patterns and other viscous flow features from the 
Navier-Stokes code solution are shown that cannot 
be easily modeled with the panel method. 

Symbols 

a speed of sound, ft/sec 

Cp pressure coefficient, 

l fuselage length, 10.328 ft 

Moc Mach number, ~ 

p pressure, lb/ft 2 

q dynamic pressure, \pV 2 , lb/ft 2 

R helicopter rotor radius, 5.164 ft 

iVR e Reynolds number, ^ 

V free-stream velocity, ft/sec 

X, Y,Z Cartesian coordinates, ft (fig. 1) 
a angle of attack 

p fluid density, slugs/ft 3 

v kinematic viscosity, ft 2 /sec 

£, t /, C curvilinear coordinates (fig. 2) 

Introduction 

Given the wide range of flight conditions in which 
helicopters must operate, particularly during hover- 
ing maneuvers, and given design constraints based 
on internal cargo and external stores, the aero- 
dynamic optimization of the fuselage is not always 
possible. However, the fuselage can significantly af- 
fect the overall performance of the helicopter in all 
flight conditions. Understanding and predicting the 
aerodynamics of helicopter fuselages will be impor- 
tant to future designs, particularly when the designs 
require greater range and speed. 

Analytical methods for evaluating the aero- 
dynamics of helicopter fuselages are available, in- 
cluding both potential theory and Navier-Stokes so- 
lutions. Early computational methods such as in 


reference 1 were based on the solution of the potential 
equation using a singularity method with constant- 
strength source panels. Since that early work, the 
computation of flow over arbitrarily shaped bodies 
has advanced significantly. Many examples of the 
extent to which panel methods have advanced can 
be found in reference 2. 

The shape of most helicopter fuselages as well as 
the wide range of flight conditions virtually guaran- 
tees that some amount of flow separation will occur. 
A computational method could model this separa- 
tion in panel methods with a boundary layer model 
(coupled inviscid and viscous solver). However, if 
separation does occur, the code must also model the 
wake. This modeling is done by shedding a wake that 
convects downstream the vorticity released when the 
boundary layer separates. The success of this ap- 
proach depends on the ability to correctly calculate 
both where the wake leaves the fuselage and its tra- 
jectory. One approach is to test the configuration 
in a wind tunnel and determine the separation loca- 
tion experimentally. This information can then be 
used in the potential code to determine the wake 
location. (See refs. 3 and 4.) More sophisticated 
approaches determine the wake separation point as 
part of the boundary layer solution. (See ref. 5.) 
Although much has been done to improve potential 
methods for computing separation, the calculation of 
helicopter fuselage flows remains challenging. 

New methods are becoming available that promise 
better predictions of complex helicopter fuselage 
flows, particularly of separation. The Navier-Stokes 
equations are the basis for computing the flow in 
complex separation regions. In the past several years, 
much progress has been made in the solution meth- 
ods needed to solve the Navier-Stokes equations. Ap- 
plications of thin-layer Navier-Stokes solvers to air- 
foils and wfngs can be found in references 6 and 7. 
In reference 8, solutions are shown for low-speed con- 
ditions over a prolate spheroid. However, few refer- 
ences are available that demonstrate the calculation 
of flow over a helicopter fuselage. Narramore and 
Brand (ref. 9) used a thin-layer Navier-Stokes code 
to study the flow over the fuselage of a Bell 214ST 
helicopter and made comparisons with experimental 
results. 

This paper uses experimental data from a generic 
helicopter fuselage shape (ref. 10) to assess these two 
methods of computation. The data used in this study 
were obtained at a Mach number of 0.062 and an 
effective Reynolds number of 4.46 x 10 6 at angles 
of attack of —10°, —5°, 0°, and 5°. By computing 
the flow over this fuselage at these angles of attack 
and by comparing both pressure distributions and 



flow features, we will summarize the strengths and 
weaknesses of these two methods for predicting this 
flow. 

Codes 

The two computational methods used for this 
study represent different approaches to the model- 
ing of fluid flow. The first method, the VSAERO 
code (ref. 11), uses potential theory with a boundary 
layer calculation coupled with the inviscid solution. 
The second method, the CFL3D code (ref. 8), solves 
the thin- layer Navier- Stokes equations. Both meth- 
ods allow the study of separation and vortical flow; 
however, the potential- theory code requires empirical 
knowledge of the separation and the Navier-Stokes 
method computes separation from first principles. 

The VSAERO code is a commercially available 
potential- theory panel-method code capable of com- 
puting flow over bodies of arbitrary shape. The body 
is represented by panels on which the source and 
doublet strengths are determined. Two-dimensional 
boundary layer calculations can be made along sur- 
face streamlines. The effects of the boundary layer 
calculations can then be coupled with the potential- 
theory solution. Wake panels can be used to simulate 
separation (e.g., behind bluff bodies or at the trailing 
edges of wings). The boundary layer calculation will 
indicate where separation will occur; however, the 
user is responsible for determining the starting loca- 
tion of the wakes. Iterations can then be performed 
on the wakes to allow them to deform to equilibrium. 
Advantages of the code include speed, ability to rep- 
resent complex geometries, and ease of use. However, 
modeling separation (e.g., behind bluff bodies) can 
be difficult. 

The CFL3D code was developed at NASA 
Langley Research Center and solves the thin-layer 
Navier-Stokes equations. The code uses a third-order 
upwind-biased method with Roe flux-difference split- 
ting to solve the equations. A multigrid scheme is 
used to improve the convergence time. The code 
is also capable of using multiblock grids, although 
for this study a single-block grid was used. Turbu- 
lence is modeled after the approaches of Baldwin and 
Lomax. (See ref. 12.) The cases presented here in- 
volve some amounts of separation where the Baldwin- 
Lomax model is uncertain (ref. 13); however, during 
these calculations this turbulence model was the only 
one available. 

Geometry 

The geometry chosen for this study was the ro- 
tor body interaction fuselage (ROBIN) that has 


been used in several helicopter investigations in the 
Langley 14- by 22-Foot Subsonic Wind Tunnel. (See 
ref. 13.) Figure 1 shows a computer simulation of the 
top and side views of this fuselage, the pressure tap 
locations used in a specific test, and the model itself 
installed in the tunnel. This body is defined analyti- 
cally so easy refinements can be made to the geome- 
try during grid development. Experimental data are 
available in the form of steady pressures at 14 sta- 
tions along the fuselage. The data given are for the 
fuselage with the rotor hub but with the blades re- 
moved. No attempt was made to model the rotor 
hub in any of the calculations. 

A C-0 volume grid was used for the configuration 
and is shown with the surface grid in figure 2. Several 
grid refinements were used to obtain the proper 
definition of the fuselage for the Navier-Stokes code. 
These refinements altered the distribution of the 
surface grid to reduce the solution dependency on the 
grid. The final volume grid consisted of 145 points 
in the streamwise direction, 65 radial points, and 
65 points in the normal direction. The surface grid 
consisted of 129 streamwise points and 65 radial 
points. 

The Gridgen code (ref. 14) was used to generate 
a single- block grid for these calculations. A single- 
block grid was chosen because future pilot Navier- 
Stokes codes may not handle multiblock grids. How- 
ever, a single-block grid makes the generation of a 
surface grid difficult, primarily at the intersection of 
the nacelle and the fuselage. The nose and tail re- 
gions of the nacelle cause disturbances in the surface 
grid that affect the gridding on the rest of the fuse- 
lage with the Gridgen algorithm used. The solution 
to this problem was to cluster the streamwise grid 
lines that make up the nacelle in a narrow region 
ahead of and behind the nacelle. (See the detail in 
fig. 2.) This change in grid spacing can cause con- 
vergence problems due to the abrupt change in grid 
cell spacing, but it allows the surface of the nacelle 
to be accurately defined. 

The surface grid taken from the volume grid was 
used to determine the paneling for the VSAERO 
code. Grid lines were removed from the surface 
grid to obtain a coarser panel geometry. The final 
representation contained 1768 panels. Even with this 
rather coarse surface grid, the VSAERO code panel 
distribution was much finer than is typically required 
to resolve the flow for this type of geometry. 

Results 

All calculations were performed at the experimen- 
tal Mach number of 0.062. Slow convergence was ev- 
ident with the Navier-Stokes code. The CFL3D code 
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was run on the Langley Cray-2 supercomputer. Sev- 
eral central processing unit hours were required to 
complete the approximately 10 000 iterations needed 
for a converged solution. Convergence was deter- 
mined from the plotted evolution of forces. In con- 
trast, the VS AERO code calculations with stream- 
line, boundary layer, and wake calculations required 
approximately 2 hours on a Silicon Graphics 320VGX 
workstation. This time depended upon the number 
of wake and boundary layer iterations required for 
the forces to converge. 

For the VSAERO code calculations, a wake was 
shed from the tail of the fuselage. Attempts to 
shed a wake from the back of the nacelle (an area 
where separation was expected to occur) did not 
give acceptable results. The location of this wake 
was determined by estimating a separation line on 
the nacelle*. The main difficulty was the tendency 
of the wake to pass through the fuselage while it 
was being relaxed. Alternatively, the wake could be 
kept rigid; however, the complex flow in this region 
could make a guess about wake geometry misleading. 
Thus, results from the region behind the nacelle are 
unlikely to compare well with the experiment or with 
the Navier-Stokes solutions. 

In reference 10, eight experimental rotor-off cases 
were presented. Four angles of attack at two dif- 
ferent free-stream flow conditions were studied. For 
this study, comparisons were made only with the low- 
speed cases, primarily to reduce the cost of the cal- 
culations. Experimental pressures at 14 stations are 
available for each case. Figures 3-6 show compar- 
isons between experiment and theory. Test condi- 
tions and pressure locations for these figures are given 
in table I. 

Pressure 

Figure 3 shows the results for the fuselage at an 
angle of attack of -10°. The pressure coefficient is 
plotted along the ordinate and the abscissa repre- 
sents the vertical coordinate of the fuselage surface. 
The individual plots represent different longitudinal 
stations along the fuselage. 

Initially, both codes compare reasonably well with 
the experiment at most stations. An examination of 
the stations where separation is expected (figs. 3(k)~ 
3(n)) reveals the advantage of the Navier-Stokes so- 
lution. At station X/R — 1.0008 (fig. 3(k)), the 
CFL3D code predicts the separation for Z/R > 0.12 
(the region behind the nacelle). Because it is a 
potential-theory code and because this region was not 
modeled with a wake, the VSAERO code calculates 
this area as a stagnation region. Thus, we expect 


Table I. Index to Test Conditions 
(a) Correlation of angle of attack to figures 3-6 


Angle of attack 

-10° 

-5° 

0° 

5° 

Figure 

3 

4 

5 

6 


(b) Correlation of pressure locations to figures 3-6 



X/R 


Part 

of figure 


0.0517 a 

0.0941 b 

0.1450 c 


0.2007 

0.2563 

0.3074 

0.3497 


d 

e 

f 

9 


0.4669 


h 


0.6003 


0.8809 ] 

1.0008 k 

1.1620 I 

1.3450 m 

1 .5298 n 


the pressure rise shown in figure 3(k). A comparison 
of the last three stations (figs. 3(l)-3(n)) shows the 
two codes are similarly accurate, except at the bot- 
tom of the fuselage. Because of the negative angle 
of attack and because of the flow disruption due to 
model installation (fig. 1), some separation might be 
expected along the bottom centerline of the fuselage. 
The experimental pressure coefficient begins to ap- 
proach 0 at this location, as do the results for the 
CFL3D code. The VSAERO code shows a pressure 
rise that is characteristic of stagnated flow. 
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Figure 4 shows the results for the fuselage at an 
angle of attack of —5°. This case shows a small differ- 
ence in predicted pressure between the two codes in 
the region ahead of the nacelle that is not seen at the 
other angles of attack. In this region, the VSAERO 
code prediction is closer to the experimental pressure 
than the CFL3D code. For stations at X/R > 1.0008 
(figs. 4(k)-4(n)), the CFL3D code calculates the 
surface static pressures more accurately. 

Figure 5 shows the results for the fuselage at an 
angle of attack of 0°. In most cases, the CFL3D code 
more accurately predicts the experimental pressures 
than the VSAERO code, especially aft of the nacelle. 
Again, this difference may relate to the fact that 
the expected separation from the nacelle was not 
modeled with the VSAERO code. Another factor 
that contributes to the discrepancy along the top of 
the fuselage is the disruption in the flow caused by 
the wake of the rotor shaft and hub (not modeled 
by either the panel or Navier- Stokes methods). This 
factor is most likely to contribute to a discrepancy at 
both the 0° and 5° angles of attack. 

Figure 6 shows the results for the fuselage at an 
angle of attack of 5°. The results for this case are also 
similar to the results for the previous cases. However, 
for this case, the CFL3D code indicates a drop in 
pressure near the top of the fuselage at stations 
X/R> 1.0008. This drop in pressure is indicative of 
vortex formation off the fuselage surface. The flow 
field solution from the CFL3D code can be visualized 
for confirmation of this vortex. 

Flow Features 

As part of the postprocessing of the Navier- Stokes 
results, particle traces were used to study the flow 
around the fuselage. Particle tracing allows the de- 
termination of separation by showing the conver- 
gence of streamlines. Particle traces, confined to 
the layer of the grid adjacent to the surface, sim- 
ulate what might be seen experimentally with oil 
flows. The contours of the normalized stagnation 
pressure were used to study flow characteristics off 
the fuselage. 

Figures 7(a)- 7(d) show the normalized stagnation- 
pressure contours as well as surface streamlines 
for the four angle- of- at tack cases calculated by the 
CFL3D code. The normalized stagnation-pressure 
contours can indicate that vortical flow is present; the 
surface streamlines will show where separation oc- 
curs. Figure 7(a) shows the relatively benign flow for 
the case at —10°. Separation is evident at the back 
of the nacelle, as expected. From the normalized 
stagnation- pressure contours, vortical flow is not in- 
dicated at the bottom rear of the fuselage. However, 


the surface streamlines indicate that a separation line 
is present along this surface. 

The case at —5° is also relatively benign, as shown 
in figure 7(b). A separation line appears along the 
lower rear portion of the fuselage that is charac- 
teristic of vortical flow. However, the normalized 
stagnation-pressure contours did not indicate a clear 
vortex that is shed from the fuselage. 

Figure 7(c) shows the case at 0°. A separation 
line can be seen along the upper rear portion of 
the fuselage. The normalized stagnation-pressure 
contours indicate that vortical flow is occurring in 
this region. 

Figure 7(d) shows the results from the CFL3D 
code for an angle of attack of 5°. On the aft portion 
of the fuselage, a vortex is shed as evidenced by the 
normalized stagnation-pressure contours. These con- 
tours also show the formation of vortical flow near the 
intersection of the nacelle and the fuselage as well as a 
separation line along the upper rear of the fuselage. 
In previous experimental tests, researchers did not 
look for these features; therefore, verification is not 
possible without further experimental measurement. 

Calculation of both on- and off-body streamlines 
is possible with the VSAERO code. The on-body 
streamlines were calculated and were used for the 
boundary layer calculation. Boundary layer calcu- 
lations are performed in a two-dimensional manner 
along these streamlines; thus, the streamlines will 
not show the flow deformation that results from the 
separation on the fuselage. Although the streamlines 
calculated by the VSAERO code will not show the ef- 
fect of separation or of vortical flow, they do provide 
a useful comparison to the Navier-Stokes results. 

Surface streamlines calculated by the VSAERO 
and CFL3D codes are shown for comparison in fig- 
ures 8(a)-8(d). The results from the VSAERO code 
do not indicate separation; however, they do allow 
some understanding of the aerodynamic character- 
istics of the fuselage. Aside from the differences 
caused by vortical flow, the main difference in the 
surface streamlines between the codes seems to oc- 
cur at locations past the nacelle. This difference is 
most likely attributable to the ability of the CFL3D 
code to calculate separation at the rear of the na- 
celle. Nevertheless, ahead of the nacelle, the two 
codes agree well. 

Conclusions 

Calculations of the flow properties over a generic 
helicopter fuselage have been presented and com- 
pared with experimental data. Potential- theory and 
Navier-Stokes methods were used for calculations at 
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four experimental conditions. Both methods agree 
well with the experiment. Prediction of flow features 
such as vortical flow and separation is highlighted. 
Relatively quick solutions are possible with poten- 
tial theory, although the ability to calculate regions 
of separation is unsatisfactory. Although more com- 
putationally expensive, the Navier- Stokes method 
allows separation and vortical flow to be studied 
from first principles. Although the Navier- Stokes 
codes predict helicopter fuselage flow sufficiently well 
for understanding flow characteristics related to vis- 
cous properties of the fluid near the body, models 
for the complex lifting rotor and its wake system 
must be incorporated for complete characterization 
of helicopter flows. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
April 4, 1994 
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(a) Location of pressure taps. 



L-78-4383 

(b) Model installed in 14- by 22-Foot Subsonic Wind Tunnel. 
Figure 1. ROBIN fuselage. 


ORIGINAL' PAGE 

BLACK AND WHITE PHOTOGRAPH 




(a) Wire frame surface detail. 



(b) C-0 volume grid. 

Figure 2. Computational grid for ROBIN. 
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Figure 3. Pressure coefficients, a = -10°; Nr- = 4.46 x 10 6 ; Afro = 0.062. 
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Figure 3. Continued. 
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Figure 3. Continued. 
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Figure 3. Concluded. 
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Figure 4. Continued. 
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Figure 4. Continued. 
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Figure 5. Continued. 
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Figure 5. Continued. 
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Figure 6. Pressure coefficients, a = 5°; Nv e = 4.46 x 10 6 ; Moo = 0.062. 
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Figure 6. Continued. 


22 



















REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


p , ^ I \sfvto i vo. u/u^-uioo 

gathering and maintaining the data ^ceTed^Jnd^ ill^aton^S^dl in , str “ ctions - searching existing data sources, 

collection of information including suggestions for reducing this burden, to Washington Headquarters Semcel ^ th,s . bur ^ en estimate or any other aspect of this 

[ Dav, Highway, Suite 1204. Arlington ^ 22202-4302, and to the Office of ManagCen^a^ 

| 3. REPORT TYPE AND DATES COVERED 

Techn ical Memorandum 

Vs. FUNDING NUMBERS - 


| 1. AGENCY USE ONLYf/.eave btank)^ 2. REPORT DATE 

June 1994 


I 4. TITLE AND SUBTITLE 

Navier-Stokes and Potential Theory Solutions for a Helicopter 
Fuselage and Comparison With Experiment 


1 6 AUTHOR(S) 

Mark S. Chaffin and John D. Berry 


WU 505-59-36-01 
PR 023-10-01-01 


1 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center Aeroflightdynamics Directorate 

Hampton, VA 23681-0001 U.S. Anny-ATCOM 

Langley Research Center 
Hampton, VA 23681-0001 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L- 17325 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM-4566 
ATCOM-TR-94-A-013 


| 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 
and 

U.S. Army Aviation and IVoop Support Command 
St. Louis, MO 63120-1798 

1 11. SUPPLEMENTARY NOTES 

Chaffin. ViGYAN, Inc., Hampton, VA; Berry: Joint Research Program Office, USAATCOM Hamnfon VA 
m »«" ‘" PPOrt,,d ' " Par *' by ‘ he Na,i °” al and Space Administration 


| 12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 


12b. DISTRIBUTION CODE 


Subject Category 02 

1 13. ABS^^^C^jMaxirmm^OO words ) 

Lo;iL! a Y>i NaV1Cr ‘ S !°. keS C ? de al , Kl a panel method code are used to predict the flow over a generic helicoDter 

£ Na ™'- St » l “ “* “'■«<»' ”e shown that cannit £ ’ W *S! 


14. SUBJECT TERMS 

Helicopter; Rotorcraft; Fuselage; Navier-Stokes; Computation; CFD; Panel method- 
Experiment ’ 

15. NUMBER OF PAGES 

28 

16. PRICE CODE 

A03 

if. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 

NSN 7540-01 280-5500 

18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 

19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 

20. LIMITATION 
OF ABSTRACT 


Prescribed by ANSI Std Z39-18 
298-102 






